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ABSTRACT 

Context. Constraining additional mixing processes and chemical composition is a central problem in stellar physics as their impact on 
stellar age determinations leads to biases in our studies of stellar evolution, galactic history and exoplanetary systems. In two previous 
papers, we have shown how seismic inversion techniques could be used to offer strong constraints on such processes by pointing 
out weaknesses in current theoretical models. The theoretical approach having been tested, we now wish to apply our technique to 
observations. In that sense, the solar analogues 16CygA and 16CygB, being amongst the best targets in the Kepler field, are probably 
the current most well observed stars to test the diagnostic potential of seismic inversions. 

Aims. We wish to use seismic indicators obtained with inversion techniques to constrain additional mixing processes in the structure 
of the components of the binary system 16Cyg. The combination of various seismic indicators will help to point out the weaknesses 
of stellar models and thus implies more constrained and accurate fundamendal parameters for these stars. 

Methods. First, we will use the latest seismic, spectroscopic and interferometric observational constraints in the litterature for this 
system to determine suitable reference models independently for both stars. We will then carry out seismic inversions of the acoustic 
radius, the mean density and of a core conditions indicator. These additional constraints will be used to improve the reference models 
for both stars. 

Results. The combination of seismic, interferometric and spectroscopic constraints allows us to obtain accurate reference models for 
both stars. However, we note that a degeneracy exists for these models. Namely, changing the diffusion coefficient or the chemical 
composition within the observational values could lead to 5% changes in mass, 3% changes in radius and up to 8% changes in 
age. We used acoustic radius and mean density inversions to further improve our reference models then carried out inversions for a 
core conditions indicator, denoted t u . Thanks to the sensitivity of this indicator to microscopic diffusion and chemical composition 
mismatches, we were able to reduce the mass dispersion to 2%, namely [0.96 M e , 1.0M o ], the radius dispersion to 1%, namely 
[1.1887? G , 1.2007? o ] and the age dispersion to 3%, namely [7.0,7.4], for 16CygA. For 16CygB, f„ offered a consistency check for the 
models but could not be used to reduce independently the age dispersion. Nonetheless, assuming consistency with the age of 16CygA 
could help to further constrain its mass and radius. We thus find that the mass of 16CygB should be between 0.93 M e and 0.96 A/ 0 
and its radius between 1.08 R e and 1.10 R a 
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1. Introduction 


In a series of previous papers (Buldgen et al. (2015) and (Buld¬ 
gen et al. (submitted)), we have analysed the theoretical aspects 
of the use of seismic inversion techniques to characterize extra 
mixing in stellar interiors. Instead of trying to determine entire 
structural profiles, as was successfully done in helioseismology 
(Basu et al. 1997, 1996; Basu & Christensen-Dalsgaard 1997) 
', we make use of multiple indicators, defined as integrated 
quantities which are sensitive to various effects in the structure. 
These indicators are ultimately new seismic constraints using all 
the available information provided by the pulsation frequencies. 


In this paper, we wish to apply our method to the binary 
system 16Cyg, observed by Kepler, for which data of unprece¬ 
dented quality is available. Moreover, this system has already 
been extensively studied, particularly since the discovery of a 
red dwarf and a jovian planet in it (see Cochran et al. 1997). 


1 Also see Christensen-Dalsgaard (2002) for an extensive review on 
helioseismology. 


Using Kepler data, this system has been further constrained by 
asteroseismic studies (Metcalfe et al. 2012; Gruberbauer et al. 
2013; Mathuret al. 2012), interferometric radii have also been 
determined (see White et al. 2013) and more recently, Verma et 
al. have determined the surface helium abundance (Verma et al. 
2014) of both stars and Davies et al. (2015) have analysed their 
rotation profiles and tested gyrochronologic relations for this 
system. 

The excellent quality of the Kepler data for these stars en¬ 
ables us to use our inversion technique to constrain their 
structure. We use the previous studies as a starting point and 
determine the stellar parameters using spectroscopic constraints 
from Ramirez et al. (2009) and Tucci Maia et al. (2014), the 
surface helium constraints from Verma et al. (2014) and the 
frequencies from the full length of the Kepler mission used 
in Davies et al. (2015) and check for consistency with the 
interferometric radius from White et al. (2013). The determi¬ 
nation of the stellar model parameters is described in Sect. 
2. We carry out a first modelling process then determine the 
acoustic radius and the mean density using the SOLA technique 
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(Pijpers & Thompson 1994) adapted to the determination of 
these integrated quantities (see Buldgen et al. 2015; Reese et al. 
2012). In Sect. 3, we recall briefly the definition and purpose 
of the indicator t u and carry out inversions of this indicator for 
both stars. We then discuss the accuracy of these results. Finally, 
in Sect. 4, we use the knowledge obtained from the inversion 
technique to provide additional and less model-dependent 
constraints on the chemical composition and microscopic 
diffusion in 16CygA. These constraints on the chemical and 
atomic diffusion properties allow us to provide accurate, yet 
of course model-dependent, ages for this system, using the 
most recent observational data. The philosophy behind our 
study matches the so-called “a la carte" asteroseismology of 
Lebreton & Goupil (2012) for HD52265, where one wishes to 
test the physics of the models and quantify the consequences of 
these changes. However, we add a substantial qualitative step 
by supplementing the classical seismic analysis with inversion 
techniques. 


2. Determination of the reference model 
parameters 

2.1. Initial fits and impact of diffusion processes 


In this section, we will describe the optimization process that 
lead to the reference models for the inversions. We carried out 
an independent seismic modelling of both stars using the fre¬ 
quency spectrum from Davies et al. (2015), which was based on 
928 days of Kepler data. A Levenberg-Marquardt algorithm was 
used to determine the optimal set of free parameters for our mod¬ 
els. We used the Cles stellar evolution code and the Lose oscilla¬ 
tion code (Scuflaire et al. 2008b,a) to build the models and cal¬ 
culate their oscillation frequencies. We used the CEFF equation 
of state (Christensen-Dalsgaard & Daeppen 1992), the OPAL 
opacities from Iglesias & Rogers (1996), supplemented at low 
temperature by the opacities of Alexander & Ferguson (1994) 
and the effects of conductivity from Potekhin et al. (1999) and 
Cassisi et al. (2007). The nuclear reaction rates we used are those 
from the NACRE project (Angulo et al. 1999) and convection 
was implemented using the classical, local mixing-length the¬ 
ory (Bohm-Vitense 1958). The empirical surface correction from 
Kjeldsen et al. (2008) was not used in this study. The cost func¬ 
tion used to perform the minimization uses the well-known for¬ 
mula: 


X 


2 


1 

N-M 


L 

i 
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where A‘ bs is an observational constraint, A‘ theo the same quan¬ 
tity generated from the theoretical model, cr, is the observational 
error bar associated with the quantity A', , N the number 
of observational constraints, and M is the number of free 
parameters used to define the model. We can already comment 
on the use of the Levenberg-Marquardt algorithm, which is 
inherently a local minimization algorithm, strongly dependant 
on the initial values. In the following section, particular care was 
taken to mitigate the local character of the results since at least 
35 models were computed independently for each star, using 
various observational constraints and initial parameter values. 
As far as the error bars are concerned, we look at the dispersion 
of the results with changes in the physical ingredients rather 
than the errors given by the Levenberg-Marquardt algorithm. 



Fig. 1: Upper panel: Observational echelle diagram of 16CygA; 
Lower panel: Observational echelle diagram of 16CygB; both 
plots illustrate the quality of the Kepler data for these stars. 
These are based on the frequencies obtained in Davies et al. 
(2015). 


We wish to emphasize that the use of other algorithms to 
select a reference model do not reduce the diagnostic potential 
of the inversions we will describe in the next sections. Indeed, 
inversions take a qualitative step beyond forward modelling 
techniques in the sense that they explore solutions outside of 
your initial model parameter space. We used various seismic 
and non-seismic constraints in our selection process and we fo¬ 
cussed our study on the importance of the chemical constraints 
for these stars. Indeed, there is a small discrepancy in the 
literature. In Verma et al. (2014), a less model-dependent glitch 
fitting techniques was used to determine the surface helium 
mass fraction, Yf. It was found to be between 0.23 and 0.25 for 
16CygA and between 0.218 and 0.26 for 16CygB (implying an 
initial helium abundance. To, between 0.28 and 0.31). In the 
seismic study of Metcalfe et al. (2012), various evolutionary 
codes and optimization processes were used and the initial 
helium abundance was 0.25 + 0.01 for a model including micro¬ 
scopic diffusion. In fact, the seismic study of Gruberbauer et al. 
(2013) already concluded that the helium mass fraction had to 
be higher than the values provided by Metcalfe et al. (2012), 
which could result from the fact that they used 3 months of 
Kepler data for their study. Therefore, the starting point of 
our analysis was to obtain a seismic model consistent with 
the surface helium constraint from Verma et al. (2014) and the 
metallicity constraint from Ramirez et al. (2009). We started by 
searching for a model without including microscopic diffusion, 
and therefore the final surface abundances Yf and Z/ are equal 
to the initial abundances To and Zo. The metallicity can be 
determined using the following equation: 


£]-*(§)-*(§).■ 


( 2 ) 


where (§) 0 is the solar value consistent with the abundances 
used in the spectroscopic differential analysis. We point out that 
in the spectroscopic study of Ramirez et al. (2009), the “solar” 
references were the asteroids Ceres and Vesta. In this study. 
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Table 1: Summary of observational properties of the system 
16CygA and 16CygB considered for this study. 



16 Cyg A 

16CygB 

R( Rg) 

1.22 + 0.02 

1.12 + 0.02 

Teff.spec (K) 

5830 + 7 

5751 +6 

Teff.phot (K) 

5839 + 42 

5809 ± 39 

Lq (Lq) 

1.56 + 0.05 

1.27 + 0.04 

[Fe/H] (dex) 

0.096 

0.051 

Yf (dex) 

[0.23,0.25] 

[0.218,0.260] 

< Av > 0 YHz ) 

103.78 

117.36 


Table 2: Optimal parameters obtained for 16CygA. 


Table 3: Optimal parameters obtained for 16CygB. 



Sb,i 

Sb ,2 

Sfl,3 

Mass (M 0 ) 

1.00839 

0.976542 

0.942557 

Radius (R 0 ) 

1.12292 

1.107 

1.09847 

Age (Gyr) 

8.16178 

7.71671 

7.37336 

T eff (K) 

5749 

5742 

5739 

Fq (Lq) 

1.23641 

1.1955 

1.17449 

Z 0 

0.0151 

0.0173 

0.0185 

To 

0.24 

0.273 

0.3015 

<+MLT 

1.567 

1.60264 

1.61543 

D 

0.0 

0.5 

1.0 

< Av > (juHz) 

117.36 

118.00 

117.37 

x 2 

0.81 

0.85 

0.88 



Sa,i 

Sa .2 

S+,3 

Mass [Mgj 

1.0523 

1.02494 

1.00146 

Radius (R 0 ) 

1.24 

1.22898 

1.21838 

Age (Gyr) 

8.23194 

7.78380 

7.33473 

Teff (K) 

5825 

5802 

5801 

Lq (Lq) 

1.58892 

1.53620 

1.50828 

Zo 

0.0165 

0.019 

0.0205 

To 

0.24 

0.271 

0.2945 

ttlVlLT 

1.61757 

1.63975 

1.67229 

D 

0.0 

0.5 

1.0 

< Av > (/;//-) 

103.74 

103.79 

103.98 

x 1 

1.18 

1.19 

1.3 


we used the (f ) 0 value from AGSS09 (Asplund et al. 2009) to 
determine the value of the metallicity Z. From the error bars 
provided on these chemical constraints, we can determine a 
two dimensional box for the final surface chemical composition 
of the model (which is the initial chemical composition if the 
model does not include any extra mixing). A summary of the 
observed properties for both components is presented in table 
1 as well as an observational echelle diagram in Fig 1 for both 
targets, illustrating the excellent quality of the seismic data for 
these stars. An initial reference model without microscopic 
diffusion was obtained using the effective temperature, T e g, the 
arithmetic average of the large frequency separation < Av > 
and the individual small frequency separations 6v„j. We did 
not include individual large frequency separations as these 
quantities are sensitive to surface effects in the frequencies and 
they would have dominated our cost-function. This would have 
been unfortunate since we want to focus our analysis on core 
regions. As we see from table 2, this model was also able to fit 
constraints such as the interferometric radius from White et al. 
(2013) and the luminosity from Metcalfe et al. (2012) although 
these quantities were not included in the x 2 of the original fit. 
The agreement between the observed and theoretical seismic 
constraints is illustrated in Fig. 2. These results might seem 
correct, but since we did not even include microscopic diffusion, 
we should consider this model as rather unrealistic in terms of 
mixing processes 2 . Therefore, we computed a few supplemen¬ 
tary models assuming a final surface chemical composition of 
Yf = 0.24 and {^) f = 0.0222 which included microscopical 

diffusion following the prescriptions of Thoul et al. (1994). In 
this case, the fit was carried out using 5 free parameters, the 
mass, the age, the mixing length parameter, denoted ctmlt, 
the initial hydrogen abundance, denoted Xq and the initial 


2 One should note that we do not imply here that microscopic diffusion 
is the only mixing process needed in a “realistic model”. 
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Fig. 2: Upper panel: Fits of the small frequency separations §02 
and §13 for 16CygA. Lower panel: Same as the upper panel for 
16CygB. (Colour online) The observational values are the green 
symbols with error bars, the red symbols are associated with 
models including solar calibrated diffusion and the blue 
symbols are associated with models without diffusion. 


metallicity, denoted Zo. We used the same constraints as for 
the first fit without diffusion, supplemented by the constraints 
on the surface chemical composition, Yf and (Z/X)f providing 
direct and strong constraints on the initial chemical composition. 

The effect of diffusion was mainly to reduce the mass, age 
and radius of the model, as illustrated in Fig. 3. This plot illus¬ 
trates the effects of diffusion for various chemical compositions 
and diffusion velocities. The subscripts 0.0, 1.0 and 0.5 are 
respectively related to a model without diffusion, with standard 
diffusion velocities and with half of these velocity values. We 
will denote this factor D in the tables presenting the results. 
Each color is associated with a particular surface chemical 
composition of these stars. All these models were fitted to 
the observed frequencies of their target. Therefore, the effect 
observed here is related to the impact of diffusion for a given 
model associated with a given set of frequencies. It is obvious 
that the reduction of the mass and radius are correlated since 
the mean density is kept nearly constant through the fit of the 
average large frequency separation. Therefore, the conclusion 
of this preliminary modelling process is that we obtain a 
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Fig. 3: Effect of the progressive inclusion of diffusion in a 
model of 16CygA. Each model still fits the observational 
constraints. 


degeneracy, meaning that we could build a whole family of 
acceptable models, inside the box of the chemical composition, 
with or without diffusion, that would be acceptable. This implies 
important uncertainties on the fundamental properties, as can 
be seen from the simple example in Fig. 3 for 16CygA. In 
the following section, we will see how the use of inversion 
techniques and especially the inversion of f„ can help us reduce 
this degeneracy and restrict our uncertainties on fundamental 
properties. Even when considering diffusion based on the 
work of Thoul et al. (1994), one should note that the diffusion 
velocities are said to be around 15 - 20% accurate for solar 
conditions. Therefore, in the particular case of 16CygA, for 
which we have strong constraints on the chemical composition, 
one can still only say that the mass has to be between 0.97 M 0 
and 1.07M Q , that the radius has to be between 1.185/? 0 and 
1.23/A; and that the age has to be between 6.8 Gy and 8.3Gy for 
this star. In other words, we have a 5% dispersion in mass, 3% 
in radius and 8% in age. 


2.2. Inversion of acoustic radii and mean densities 

In this section, we will shortly present our results for the 
inversion of the mean density and the acoustic radius. First, 
we note that the inverted results for the mean density and the 
acoustic radius are slightly different from each other. There 
is a dispersion of around 0.5% for both p and r depending 
on the reference model used for the inversion. We therefore 
consider that the results are respectively ta = 4593 ± 15 s and 
Pa = 0.83 ± 0.005 g/cm 3 to be consistent with the dispersion 
we observe. For 16CygB, we obtain similar results, namely 
tb = 4066 + 15s and f>n = 1.066 + 0.005 g/cm 3 . The kernels are 
well fitted, as can be seen for a particular example in Fig. 4. 

This justifies the fact that linear inversions are said to be 
“nearly model-independent". We emphasize that the physical 
ingredients for each model were different and that the dispersion 
of the results is smaller than 0.50%. Before the inversion, the 
dispersion of the mean density was of about 0.95% and signif¬ 


icantly different from the inversion results. In that sense, the 
model-dependency of these methods is rather small. However, 
the error bars determined by the simple amplification of the ob¬ 
servational errors is much smaller than the model-dependency, 
thus one has to consider that the result is accurate within the 
dispersion due to the reference models rather than using the 
error bars given by the inversion. Nevertheless, this dispersion is 
small and therefore these determinations are extremely accurate. 

We also observed that including additional individual large 
frequency separations in the seismic constraints could improve 
the determination of both the acoustic radius and the mean 
density of the model. However, this can reduce the weight 
given to other seismic constraints and as we will see in the next 
section, we can improve the determination of reference models 
using directly the acoustic radius and the mean density as 
constraints in the fit. We also note that neither the mean density, 
nor the acoustic radius could help us disentangle the degeneracy 
observed in the previous section for the chemical composition 
and the effects of diffusion. Indeed, these quantities are more 
sensitive to changes in the mixing-length parameter, (Xmlt, or 
strong changes in metallicity. However, as will be described in 
the following section, they can be used alongside other inverted 
structural quantities to analyse the convective boundaries and 
upper layers of these stars. 


2.3. Determination of new reference models 

After having carried out a first set of inversions using the acous¬ 
tic radius and the mean density, we carried out a supplementary 
step of model parameter determination, replacing the average 
large frequency separation by the acoustic radius and the mean 
density themselves. We obtained a new family of reference 
models which were slightly different from those obtained using 
the average large frequency separation. 

If we compare the model parameters obtained using r and 
p for the model with Yf = 0.24 and {Z/X)f = 0.0222 with 
those obtained with < Av >, presented in table 2, we note that 
there is a tendency to slightly reduce the mass and to increase 
the mixing length parameter. What is more surprising is that 
when computing individual frequency differences between the 
observed stars and the reference models, we see that using the 
acoustic radius and the mean density allows us to obtain signifi¬ 
cantly better individual frequencies. This is a by-product of the 
use of inversion techniques that could be used to characterize 
stars in a pipeline such as what will be developed for the coming 
PLATO mission. 

Considering that these models are improved compared to 
what was obtained using the large frequency separation, we 
compute a family of models for different values of Yf and (f ) / .- 

For each particular chemical composition, we computed models 
with and without microscopic diffusion. The properties of some 
models of this family are summarized in table A.l. As can be 
seen, some of the models do not reproduce well the results for 
the effective temperature or the interferometric radius. This 
means that we can use non-seismic constraint as indicators of 
inconsistent models in our study, although one should be careful 
on the conclusions derived from these quantities. For instance, 
the interferometric radii are different from the radii computed 
with the Cles models and some differences might result from 
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Position r/R 



0 0.2 0.4 0.6 0.8 1 
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Fig. 4: (Colour online) Upper panel: Example of Kernel fits for the inversion of the acoustic radius of 16CygA (Averaging kernel 
on the left and cross-term kernel on the right). Lower panel: Kernel fits for the inversion of the mean density for 16CygA 
(Averaging kernel on the left and cross-term kernel on the right). The target functions are in green and the SOLA kernels in blue. 


the very definition of the radius. One should also note that these 
results are not totally incompatible since White et al. (2013) 
conclude that the radius of 16CygA is 1.22 ± 0.02 R Q and we 
find values around 1.185 and 1.23, outside the lcr errors for 
the lower part of our dispersion. The stellar luminosity is also 
dependent on these radii values and therefore should be consid¬ 
ered with care. At the end of the day, the effective temperature 
can be constraining although there might be a slight difference 
stemming from discrepancies between the physical ingredients 
in the stellar atmosphere models used for the spectroscopic 
study of Ramirez et al. (2009) and Tucci Maia et al. (2014) and 
those used in the Cles models in this paper. However, the incon¬ 
sistencies observed for certain of these models are too important 
and therefore they should be rejected. The combination of all 
the information available will be described in Sect. 4. In the 
next section, we will use these models as references for our 
inversions of the t u indicator. One should note that this first step 
was beneficial since obtaining reference models as accurate as 
possible for these stars is the best way to obtain accurate results 
for the more difficult inversion of the f„ indicator. 


3. Inversion results for the t u core conditions 
indivator 

3.1. Definition of the indicator and link to mixing 
processes 

In Buldgen et al. (submitted), we defined and tested a new indi¬ 
cator for core conditions, applicable to a large number of stars 3 
and very sensitive to microscopic diffusion or chemical compo¬ 
sition mismatches in the core regions between the target and the 
reference model. The definition of this quantity was the follow¬ 
ing: 



where u is the squared isothermal sound speed, defined as u = £ 
f(r ) is a weighting function defined as follows: 


f(r) — r(r — R) 2 exp 7 ^ j j . (4) 

Due to the effects of the radius differences between the observed 
target and reference model, we noted that the quantity measured 
was ^f-, where R„ lr is the target radius. In Fig. 5, we illustrate 
the changes of the quantity due to the effects of diffusion for 
two of our reference models, having the same surface chemical 
composition and fitting the same observational constraints. One 
can also see the effect of surface helium and metallicity changes 

3 Provided that there is sufficient seismic information for the studied 
stars. 
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Fig. 5: Left panel: Effect of diffusion, metallicity changes and helium abundance changes on the core regions for models Sa,ci> 
Sa,c2, Sa./.i, Sa,£/i on the target function of r„, since the quantity is integrated, the sensitivity is greatly improved. Right panel: the 
Y(x ) profile of these models is illustrated, thus showing the link between t u and chemical composition and thus, its diagnostic 

potential. 


on the profile of the integrant of Eq. 3. The whole parameter set 
of these models is given in table A. 1 as well as the explanation 
of the naming convention. The diagnostic potential of the t u in¬ 
version is therefore clear, although the weighting function could 
be adapted to suit other needs if necessary. The inversion of this 
integrated quantity can be made using both the (uo, T 1 ) or the 
(uq, Y) kernels. 


3.2. Inversion results for 16 CygA 

The inversion results are summarized in Fig. 8 (they are repre¬ 
sented as orange x in the p - plot) and illustrated through an 
example of kernel fits in Fig. 6. We tried using both the (mo,Ei) 
and the (uq, Y) kernels. The high amplitude of the T 1 cross-term 
leads us to present instead the results from the (uo, Y) kernels 
although they are quite similar in terms of the inverted values. 
However, one should note that the error bars are quite important, 
and we have to be careful when interpreting the inversion results. 

This effect is due to both the very high amplitude of the 
inversion coefficients and the amplitude of the observational 
error bars. When compared to the somewhat underestimated 
error-bars of the acoustic radius and mean density inversion, it 
illustrates perfectly well why it is always said that two inversion 
problems can be completely different. In this particular case, 
using various reference models allows us to already see a trend 
in the inversion results. We clearly see that the value of f„ for 
our reference models is too low and we see that the dispersion 
of the inversion results is rather low, despite the large error bars. 
One should also note that the quality of the kernel fit is also a 
great indicator of the quality of the inverted result. For most 
cases, the kernels were very well fitted and the low dispersion of 
the results means that there is indeed information to be extracted 
from the inversion. We will see how this behaviour is clearly 
different for 16CygB. 

Nevertheless, one could argue that a small change of t u 
could be easily obtained through the use of diffusion, or 
chemical composition changes. We will see in Sect. 4 how 
combining all the information along with new constraints from 


the inversion technique can be extremely restrictive in terms of 
chemical composition and diffusion processes. Indeed, t u should 
not be considered as a model-independent age determination 
or as an observed quantity disentangling all physical processes 
occurring in stellar cores. In fact, it is simply a nearly model- 
independent determination of a structural quantity optimized to 
more sensitive to any change in the physical conditions in stellar 
cores than classical seismic indicators. The amplitude of the 
error bars remind us that this sensitivity comes at a cost and in 
this study we will consider that having a reference model with 
a -k- ~ 3.2 or 3.3 will be acceptable if it still fits the other 

observational constraints. 


3.3. Inversion results for 16 CygB 


The case of 16CygB is completely different. In fact, whereas 
the inversion for the acoustic radius and the mean density have 
been successful and we could build improved models for this 
star, the inversion of the t u indicator was less successful. The 
results were good, in the sense that the kernels are well fitted. 
However, we can see from Fig. 7 that the amplification of the 
observational error was too high to constrain the microscopic 
diffusion effects or the chemical composition. In fact, it is not 
surprising since the error bars on the observed frequencies are 
larger than for 16CygA. 

As a matter of fact, the observational error dominate the 
inversion result, as can be easily shown in Fig. 7. We see that the 
relative change in t u is smaller when microscopic diffusion is 
included in the model but this is due to the fact that the inversion 
result is closer to the reference value rather than the opposite. 
Thus, this means that t u can be used as a consistency check for 
future investigations, to ensure that we stay within the error 
bars of the inverted value, but we cannot hope to gain additional 
information for this star from this indicator. 
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Fig. 6: (Colour online) Example of kernel fits for the t u inversion. The left panel is associated with the averaging kernel and the 
right panel is associated with the cross-term kernel. The target functions are in green and the SOLA kernels in blue. 
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Fig. 8: (Colour online) Results of the t u inversions for 16CygA (left panel) and positions of the reference models in the chemical 
composition box derived from spectroscopic and seismic constraints (right panel). The orange x are the inversion results whereas 
the other symbols are associated with various reference models the positions of which are shown in the right-hand side plot. The 
colour is associated with the Yf value, the type of symbol with the value, and the size of the symbol with the inclusion of 

diffusion. 


4. Constraints on microscopic diffusion and 
chemical composition 

4. 1. Reducing the age , mass and radius dispersion of 
16 Cyg A 

In this section, we will use the information given by t u to further 
constrain chemical composition and microscopic diffusion. 


Previously, we always ensured that the reference models were 
inside the chemical composition box that was defined by the 
constraints on surface helium obtained by Vermaetal. (2014) 
and the spectroscopic constraints obtained by Ramirez et al. 
(2009). In Sect.3.2, we concluded that our model should have 
at least a ~ 3.2 or 3.3 or higher. The first question that 

arises is whether it is possible to obtain such values for ^ given 
the constraints on chemical composition. The second question 
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Fig. 7: 4 inversion results for 16CygB. (Colour online) The red 
+ are the reference models and the blue x the inverted results. 
The lower + are associated with the upper x and refer to models 
including solar-calibrated diffusion. 


is related to the impact of microscopic diffusion. 

In fact, t u is a measure of the intensity of the uq gradients 
in the core regions. Thus, since uq ~ j r where T is the tem¬ 
perature and p the mean molecular weight, including diffusion 
will increase the p gradients, since it leads to the separation 
of heavy elements from lighter elements. It is then possible 
to increase significantly the diffusion speed of the chemical 
elements and to obtain a very high value of f„ for nearly any 
chemical composition. However, in Thoul et al. (1994), the 
diffusion speed are said to be a; 15 - 20% accurate, and suited 
for solar conditions. Moreover, since increasing diffusion also 
accelerates the evolution, we could also end up with models 
that are too evolved to fit simultaneously t u , the chemical 
composition constraints and the seismic constraints. Looking 
at the parameters of our reference models, we note that we 
are indeed very close to solar conditions and we suppose that 
our diffusion speed should not be amplified or damped by 
more than 20%. The results of this analysis are summarized in 
Fig. 8, which is a p - plot where the reference models and 
the inverted results are represented. In what follows, we will 
describe more precisely our reasoning and refer to Fig. 8 when 
necessary. We used a particular colour code and type of symbol 
to describe the changes we applied to our models. Firstly, colour 
is associated with the final surface helium mass fraction Yf, 
namely, blue for Yf = 0.24, red if Yf < 0.24 and green if 
Yf > 0.24. Secondly, the symbol itself is related to the (§)^ 

namely, a x for (j) < 0.0222, a o for ( jf'j = 0.0222 and a o for 

(f) > 0.0222. The size of the symbol is related to the inclusion 

of microscopic diffusion, for example the large blue and red dots 
in Fig. 8 are related to models including microscopic diffusion. 

Since increasing diffusion should increase the t u value, we 
compute a model, with Yf - 0.24 and = 0.0222, including 

diffusion from Thoul et al. (1994), fitting the seismic constraints 
and the effective temperature. This model is represented by the 


Table 4: Accepted parameters obtained for 16CygA when 
taking into account the constraints from the inversion of t u . 



Accepted 16CygA models 

M (Mg) 

0.96-1.00 

Age (Gy) 

7.0-7.4 

To (dex) 

0.30-0.31 

Z 0 (dex) 

0.0194-0.0199 

D 

1.00-1.15 

Q'mlt 

1.75-1.90 

L( L 0 ) 

1.49-1.56 

R( Re) 

1.19-1.20 


large blue dot and we note that including diffusion improves 
the agreement, but is not sufficient to reach what we defined to 
be our acceptable values for This is illustrated by the fact 
that in Fig. 8, the large blue dot is above the small blue dot. 
Therefore we decide to analyse how f„ depends on the chemical 
composition. To do so, we compute a model for each corner 
and each side of the chemical composition box. These models 
are represented in Fig. 8 by the o, o and x of various colours. 
From these results, we see that increasing the helium content, 
namely considering that Yf e [0.24,0.25] increases t u , as does 
considering (j)j. e [0.0209,0.0222], In simpler terms, we see 

that the green dot and the blue cross are above the blue dot in 
Fig. 8. The first tendency is quickly understood since increasing 
the helium abundance leads to higher central p and therefore a 

local minimum in the t/o profile. Because t u is based on (jij , 
this does not imply a reduction of the value of the indicator, but 
an increase due to a secondary lobe developing exactly in the 
same way as what happens when including diffusion (see Fig. 
3). The second tendency can be understood by looking at the 
central hydrogen abundance. In this case, we see that the central 
hydrogen abundance is reduced and thus the mean molecular 
weight is increased and lead to a minimum in u tt in the centre. 
One should note that this effect is not as intense as the change in 
helium but is still non-negligible. 

Therefore, our seismic analysis favours models that lie 
within Y f € [0.24,0.25] and (§) 6 [0.0209,0.0222], Including 

diffusion in these models increases further the value and 
brings it in the range of the 3.2, 3.3 values, which is much 
more consistent with the inversion results. These final models 
are represented in Fig. 8 by the large green +. One should also 
note that an upper boundary can be drawn from the effective 
temperature and the seismic constraints. In other words, the fit 
of the other quantities can slightly increase the x~ up to values 
of 1.6 and thus slightly reduce the quality of the fit. This is not 
alarming but still means that one should not put all the weight 
of the fit of the model on the inversion results but try to find a 
compromise between seismic, spectroscopic, and inverted con¬ 
straints. Considering models that fulfil all these constraints, we 
are able to reduce the degeneracy previously observed. We thus 
conclude that the mass of 16CygA must be between 0.97M 0 
and 1.0M o and its age must be between l.OGy and 1 AGy. These 
values are subject to the hypotheses of this study and depend 
on the physics used in the stellar models (opacities, nuclear 
reaction rates, abundances). We recall here that there is no way 
to provide a seismic fully model-independent age, but inversions 
allow us at least to check the consistency of our models with 
less-model dependent structural quantities. These consistency 
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checks can lead to a refinement of the model parameters and in 
this particular case to constraints on microscopic diffusion. 

For the sake of completion, we also analysed the impor¬ 
tance of the abundances used to build the model. In fact, the 
[Fe/H] constraint being extremely dependent on the solar 
, we wanted to ask the question of whether the inversion 
would have also provided a diagnostic if we had used the 
GN93 abundances to determine the metallicity. Using these 
abundances and the associated (f ) 0 which is equal to 0.0244, 
one ends up with models having much higher metallicities, of 
the order of 0.0305 when no diffusion is included in the model. 
In fact we ended up with the same tendencies in the chemical 
composition box, but with completely different values of 
implying slightly higher masses of around 1.03 M e and slightly 
lower ages around 6.8Gy. However, when carrying out the t u 
inversion, we noted that we still had to increase the helium 
content, include diffusion, and reduce the (j^). The interesting 

point was that even the lowest (jj, associated with the highest 
Yf with increased diffusion could not produce a sufficiently 
high value of t u . In that sense, it tends to prove what we already 
suspected, that the GN93 abundances should not be used in 
the spectroscopic determination of the (jj for this study. This 
emphasizes the importance of consistency with the differential 
spectroscopic study that determined the [Fe/H] value for the 
star one wishes to study. In this particular case, we see that the 
inversion of t u is able to detect such inconsistencies, thanks 
to its sensitivity on metallicity mismatches. However, if the 
model is build with the Z/X determined from the AGSS09 
solar reference value, but using the GN93 solar heavy element 
mixture, we cannot detect inconsistencies. In fact, we obtain 
the same conclusion as before since these models are nearly 
identical in terms of internal structure. 


4.2. Impact on the mass and radius dispersion of 
16 Cyg B 

In the previous section, we used the t„ inversion to reduce the 
age, mass and radius dispersion of 16CygA. Moreover, we know 
from Sect. 3.3 that the inversion of t u for 16CygB can only 
be used to check the consistency of the model but not to gain 
additional information. However, since these stars are binaries, 
we can say that the age values of the models 16CygB must be 
compatible with those obtained for 16CygA. From the inversion 
results of 16CygA, we have also deduced that we had to include 
atomic diffusion in the stellar models and since both stars are 
very much alike, there is no reason to discard microscopic 
diffusion from the models of the B component when we know 
that it has to be included in the models for the A component. 

Therefore, we can ask the question of what would the 
mass and radius of 16CygB be if one includes diffusion as 
in 16CygA and ensures that the age of the models remain 
compatible. The question of the chemical composition is also 
important since Ramirez et al. (2009) found a somewhat lower 
value for the [Fe/H] of the B component and Verma et al. (2014) 
found larger uncertainties for the surface helium abundance, 
although the centroid value was the same as that of 16CygA. To 
build these new models, we imposed that they included atomic 
diffusion with a coefficient D of 1.0 or 1.15. The age was to be 
between 7.0 Gy and 7.4 Gy. The metallicity was imposed to 


Table 5: Accepted parameters obtained for 16CygB when 
taking into account the constraints on 16CygA. 



Accepted 16CygB models 

M(M 0 ) 

0.93 - 0.96 

Age (Gy) 

7.0 - 7.4 

To (dex) 

0.30-0.31 

Z Q (dex) 

0.0151-0.0186 

D 

1.00-1.15 

Q'mlt 

1.65-1.80 

L (L 0 ) 

1.17-1.24 

R GU) 

1.08-1.10 


be within the error bars provided by Ramirez et al. (2009) and 
the surface helium abundance was to be within [0.24,0.25]. We 
used the same constraints as before to carry out the fits using 
the Levenberg-Marquardt algorithm and found that the mass 
was to be within 0.93 M 0 and 0.96 M 0 , thus a 1.5% dispersion 
and the radius was to be within 1.08 R Q and 1.10 R Q , thus a 
1% dispersion. We would like to emphasize here that these 
values are of course dependent on the results of the modelling 
of 16CygA and are thus more model-dependent since they do 
not result from constraints obtained through seismic inversions 
but are a consequence of the binarity of the system. It is clear 
that a changes in the values of the fundamental parameters for 
16CygA will induce a change in the values of 16CygB. 


4.3. Discussion 

The starting point of this study was the determination of funda¬ 
mental parameters for both 16CygA and 16CygB using seismic, 
spectroscopic and interferometric constraints. However, the 
differences between our results and those from Metcalfe et al. 
(2012) raise questions. One could argue that the inversion leads 
to problematic results and that the diagnostic would have been 
different if the surface helium determination from Verma et al. 
(2014) would have not been available. 

Therefore, for the sake of comparison, we asked the ques¬ 
tion of what would have been the results of this study if we had 
not included the surface helium abundance from Verma et al. 
(2014) in the model selection process. We carried out a few 
supplementary fits, using the mass, age, o:mlt, Xq and Zq as 
free parameters, using all the previous observational constraints 
as well as the prescription for microscopic diffusion from 
Thoul et al. (1994), but excluding the Yf value. The results 
speak for themselves since we end up with a model having a 
mass of l.O9M 0 and an age of 7.19Gv compatible with the 
results from Metcalfe et al. (2012). This means that the deter¬ 
mining property that leads to the changes in the fundamental 
parameters of the star was, as previously guessed, the surface 
helium value. Hence, without this Yf constraint, one would end 
up with two solutions with completely different masses and 
ages but fitting the same observational constraints. This does 
not mean that the results from Metcalfe et al. (2012) are wrong, 
they simply were the best results one could obtain without the 
surface helium constraint and with 3 months of Kepler data. In 
fact, this is only an illustration of the importance of chemical 
composition constraints in stellar physics. The To - XI trend has 
already been described in Baudin et al. (2012) and the fact that 
we find lower masses when increasing the helium abundance is, 
ultimately, no surprise. 
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At this point, we wanted to know what the inversion re¬ 
sults would have been if we had used reference models with 
similar parameters as what was obtained in Metcalfe et al. 
(2012). We ended up with similar results for both the acoustic 
radius and the mean density inversion, but more interestingly, 
the t u inversion also provided non-negligible corrections for 
this model. In fact, even with microscopic diffusion, the '-A^- 

K Ref 

value was of: 2.12g/cm 6 whereas the inverted result was: 
- 3.5 + 0.5 g/cm 6 . Therefore the diagnostic potential of the 

obs * 

indicator is still clear, since it could have provided indications 
for a change in the core structure of the model. Assuming that 
diffusion velocities are around 20% accurate, one could have 
invoked either an extra-mixing process or a change in the initial 
helium composition to explain this result. Disentangling be¬ 
tween both cases would then have probably required additional 
indicators. 


5. Conclusion 

In this article, we have applied the inversion techniques pre¬ 
sented in a series of previous papers to the binary system 
16CygA and 16CygB. The first part of this study consisted 
in determining suitable reference models for our inversion 
techniques. This was done using a Levenberg-Marquardt 
algorithm and all the seismic, spectroscopic and interferometric 
observational constraints available. We used the oscillation 
frequencies from Davies et al. (2015), the interferometric radii 
from White et al. (2013), the spectroscopic constraints from 
Ramirez et al. (2009) and the surface helium constraints from 
Verma et al. (2014). 

Due to these constraints on the surface chemical composi¬ 
tion, our results are different from those of Metcalfe et al. 
(2012). The test case we made without using the constraint 
on surface helium from Verma et al. (2014) demonstrates the 
importance of constraints on the chemical composition for 
seismic studies. In fact, having to change the initial helium 
abundance from 0.25 to values around 0.3 is of course not 
negligible. This emphasizes that we have to be careful when 
using free parameters for the stellar chemical composition in 
seismic modelling. The same can be said for the constraints 
on the stellar [Fe/H] from the study of Ramirez et al. (2009). 
For this particular constraint, we have to add the importance of 
the solar mixture used in the spectroscopic study. Due to the 
important changes in the (f ) 0 from the GN93 abundances to the 
AGSS09 abundances, consistency with the spectroscopic study 
has to be ensured. Otherwise, the final results of the seismic 
modelling could be biased. In this case, we chose to use the 
AGSS09 solar mixture. We note that our reference models tend 
to be consistent with the spectroscopic, seismic and interfer¬ 
ometric constraints and that independent modelling of both 
stars leads to consistent ages. We also note the presence of a 
certain modelling degeneracy in terms of chemical composition 
and microscopic diffusion. Accordingly, we could obtain rather 
different values for the mass, the radius and the age of both stars 
by assuming more intense diffusion and changing the chemical 
composition within the error bars from both Ramirez et al. 
(2009) and Verma et al. (2014). We also note that when not 
considering the constraints on surface helium, we obtained 
results compatible with Metcalfe et al. (2012) but the t u values 


were too low even when diffusion was included in the models. 
This reinforces the importance of constraints on the chemical 
composition and illustrates to what extent inversions could be 
used given their intrinsic limitations. 

Having obtained suitable reference models, we then car¬ 
ried out inversions for the mean density , p, the acoustic radius, 
r, and a core condition indicator, t u . The first two quantities 
were used to improve the quality of the reference models. As 
a by product, we noted that models fitting both p and r were 
in better agreement in terms of individual frequencies. We also 
found that both of these quantities could not disentangle the 
effect of the degeneracy in terms of diffusion and chemical 
composition. However, they could be well suited to analyse 
uppers layers along with other quantities. 

After the second modelling process, we carried out inver¬ 
sion for the t„ indicator and noted that the degeneracy in terms 
of chemical composition and diffusion could be reduced for 
16CygA. In fact, to agree with the inverted result, one has to 
consider diffusion speed calibrated for the sun or slightly higher 
(by 10% or 15%). Values higher than 20% were considered not 
physical by Thoul et al. (1994) and were therefore not analysed 
in this study. At the end of day, we come up with a smaller 
dispersion in terms of mass and age for 16CygA, namely that 
this component should have a mass between 0.97M 0 and 1.0M 0 , 
a radius between 1.188R 0 and 1.200/? o and an age between 
7.0 Gy and 7.4Gv. Again the slight differences between the 
seismic radius provided here and the interferometric radius 
might stem from different definitions of the interferometric 
radius and the seismic one. We also conclude that the t u 
inversion for 16CygB could only be used as a consistency 
check but could not help reduce the dispersion in age. However, 
as these stars are binaries, a reduced age dispersion for one 
component means that the second has to be consistent with this 
smaller dispersion. Therefore, we were able to deduce a smaller 
mass and radius dispersion for the second component, namely 
between 0.93 M 0 and 0.96 M 0 and between 1.08 R Q and 1.10 R Q . 

Finally, we draw the attention of the reader to the following 
points: the age dispersion we obtain is not model-independent, 
we assumed physical properties for the models and assumed 
that the agreement in t u was to be improved by varying the 
chemical composition within the observational constraints 
and by calibrating microscopic diffusion. This does not mean 
that no other mixing process has taken place during the evo¬ 
lutionary sequence that could somehow slightly bias our age 
determination. In that sense, further improved studies will be 
carried on, using additional structural quantities, more efficient 
global minimization tools for the selection of the reference 
models and possibly improved physical ingredients for the 
models. In conclusion, we show in this study that inversions 
are indeed capable of improving our use of seismic information 
and therefore through synergies with stellar modellers, helping 
us build new generations of more physically accurate stellar 
models. 
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Sb,ci 

Ss,C2 

M( M 0 ) 

1.00793 

0.960578 

R (Re) 

1.10560 

1.08821 

Age (Gyr) 

8.16170 

7.23565 

Teff (K) 

5793 

5829 

Lq (Lq) 

1.23498 

1.22672 

Zo 

0.0151 

0.01805 

To 

0.24 

0.29183 

QAilt 

1.667 

1.78034 

D 

0.0 

1.0 


Appendix A: Intermediate results of the forward 
modelling process 

After the first step of forward modelling, we carried out supple¬ 
mentary fits to obtain new reference models for both 16CygA 
and 16CygB. In fact, we replaced the average large frequency 
separation by the acoustic radius and the mean density, as dis¬ 
cussed in Sect. 2.3. We used the following naming convention 
for these models: the first letter, A or B is associated with the 
star, namely 16 Cyg A or 16 CygB\ the second letter is associated 
with the chemical composition box in the right pannel of Fig. 8, 
C being the central chemical composition, L the left-hand side, 
R the right-hand side, U the upper side and D the lower side 
(D for down); the number 1 or 2 is associated with diffusion, 1 
is for models without microscopic diffusion, 2 is for models in¬ 
cluding the prescriptions of Thoul et al. (1994) for microscopic 
diffusion. These results are illustrated in the following tables for 
both stars: 
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Table A.l: Optimal parameters obtained for 16CygA using the acoustic radius and the mean density rather than < Av >. 



Sa.ci 

Sa,C2 

Sa.c/i 

Sa,;/2 

Sa.di 

Sa,D2 

Sa.ri 

Sa„R2 

Sa.li 

Sa,L2 

M(M e ) 

1.049 

0.999 



1.06 


1.055 

1.001 

1.049 

0.983 

R( Re) 

1.221 

1.201 

1.216 

1.198 

1.227 


1.222 

1.201 

1.220 

1.195 

Age (Gyr) 

8.30 

7.38 


6.77 

8.33 

7.53 

8.34 

7.31 

8.11 

7.33 

Teff (K) 

5852 

5828 


5992 

5842 

5811 

5827 

5837 

5912 

5877 

(Lo) 

1.570 

1.494 

1.61255 

1.66196 

1.574 

1.482 

1.546 

1.504 

1.633 

1.529 

Zo 

0.0165 

0.0205 


0.0195 

0.0167 


0.0174 

0.0210 

0.0155 


Y 0 

0.24 

0.2945 

IB 

0.307729 

0.23 


0.24 

0.2968 

0.24 


&MLT 

1.68 

1.74 

1.75 

1.97 

1.68957 

1.72 

1.67 

1.76 

1.75 

1.78 

D 

0.0 

1.0 


1.0 

0.0 


0.0 

1.0 

0.0 
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